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Abstract 

We report the observation of the twisted decay of quadruply charged vortices in an atomic Bose- 
Einstein condensate. Supporting numerical simulations show that the singly-charged vortices, 
which result from the decay of a multi-charged vortex, twist around intertwined in the shape of 
helical Kelvin waves. Finally, we propose to apply this effect to generate an almost isotropic state 
of turbulence which we characterize in terms of the velocity statistics. 
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Superfluids are noteworthy because they flow without dissipating energy. Even more re¬ 
markably, superfluid flow patterns are characterized by the quantization of vorticity, arising 
from the existence and the uniqueness of a macroscopic wave function 'ip. The flow velocity 
is proportional to the the phase gradient of and the circulation around a vortex line must 
be an integer multiple n = 1,2, ••• of the quantum of circulation k, = h/m, where h is 
Planck’s constant and m is the atomic mass. The singular nature of this quantized vorticity, 
concentrated along lines, has yet another important consequence: within a turbulent tangle 
of quantum vortices, the velocity components obey power-law statistics ina, unlike Gaus¬ 
sian statistics typical of ordinary turbulence. Recent experiments and numerical simulations 
j3] have shown that, under certain conditions, the turbulent superflows share a remarkable 
property with classical turbulence: the same Kolmogorov energy spectrum [1], describing 
the distribution of kinetic energy over the large length scales. This property suggests that 
quantum turbulence may be the ‘skeleton’ of classical turbulence [H]. 

Atomic Bose-Einstein condensates (BECs) are emerging as ideal systems to explore the 
quantization of vorticity and many other fundamental problems concerning the nature of 
turbulence [Hj. Vortices are more easily nucleated, manipulated [ZlIH! and observed [MI] 
in BECs than in superfluid helium, due to the typical vortex core sizes which are orders of 
magnitude larger in gaseous BECs (~ 10“^ m) than in liquid He (~ 10“^° m). 

Multiply quantized vortices are interesting especially since they are energetically unstable 
and decay into singly quantized vortices [I2HIS]. The angular momentum and the energy 
of an isolated vortex in a homogeneous superfluid grow respectively with n and nP HE]. 
Eor the same angular momentum, multi-charged (n > 1) vortices carry more energy, and, 
in the presence of dissipative mechanisms (e.g. thermal excitations), tend to decay into 
singly-charged vortices, minimizing the system’s energy. Besides the energy instability, there 
is also a dynamical instability dlHISI, which is able to destabilize the vortices even in 
the absence of dissipation, at zero temperature, whose time-scales have been investigated 
[HI (HI dSl [IHl EDH23I. The technique of topological phase imprinting [21| has allowed the 
controlled generation of multi-charged vortices [HI [II E21 EH] in atomic condensates. The 
splitting of doubly quantized vortices into two singly quantized vortices has been observed 
in a Na BEC and the decay process was studied HaEOlEl]. Quadruply quantized vortices 
are theoretically predicted |26| to decay presenting various symmetries of splitting, making 
them an interesting research topic. Recent work has determined that the stability of such 
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vortices is affected by the condensate’s density na and size EH. and by the nature of the 
perturbations [22! • 

In this letter, we discuss the decay of an n = 4 charged vortex line nucleated in pro¬ 
late ®^Rb BECs. By combining experimental observations with numerical simulations, we 
present results showing the twisted decay process of a multi-charged vortex, induces helical 
Kelvin waves on the resulting singly-charged vortex lines. We also show that the observed 
intertwined decay of multi-charged vortices may be exploited to create an almost isotropic 
state of quantum turbulence in atomic condensates. Kelvin waves are of particular interest 
because they are thought to play a key role in quantum turbulence 1271, being observed in 
classical fluids and superfluid helium, including atomic condensates [281 [29| . 

The experimental sequence to produce the BEC runs as follows. First, ®^Rb BECs are 
produced in the \F = 2, mp = -\-2) hyperfine state, with a small thermal fraction in a cigar¬ 
shaped QUIC magnetic trap. The atoms undergo forced evaporation for 23 s, from 20 MHz to 
about 1.71 MHz, following a non analytic curve, experimentally determined by optimizing 
the phase space density and the runaway condition in each step. We typically produce 
samples with no more than 35% of thermal fraction, well into the Thomas-Fermi regime. 
The measured trapping frequencies are: = 2ttx 21.1(1)Hz in the symmetry axis, and u)r = 

27r X 188.2(3)Hz in the radial direction, resulting in the (geometric) mean trap frequency, 
(h = 27r X 90.7(4)Hz, and the harmonic oscillator length, Uho = \/h/muj ~ 1.13/rm. The 
typical properties of our BECs are: reduced temperature T/Tq ss 0.65—0.70; A^o ~ 1—2 x 10^ 
condensate atoms; ~ 1 — 3 x 10® total number of atoms; healing length ^ ~ 0.17/im; 
axial and radial Thomas-Fermi radii Rtf( 2:,0) = 13.9(7)/rm and = 5.8(6)/rm, 

respectively; chemical potential /r ~ 100 nK; anisotropy parameter A = oJz/oJr ~ 0.112. The 
multi-charged vortices are nucleated by adiabatically reducing and inverting the initial trap 
bias field from B^it = 0) ~ 0.5G down to = r) ~ —0.5G, in r ~ 5ms typically. 

The resulting magnetic bias field along the weak trap direction, is reversed during the 
process [HD]- We find that the n = 4 charged vortices, topologically imprinted, are unstable 
and decay into four singly-charged vortices in a twisted unwinding manner. The process 
is very reproducible, over several runs and different time-of-ffight values. The 2D optical 
depth (OD) images presented correspond to the integrated number density acquired along 
the symmetry axis (axial images), shown in the top of Fig. or integrated along y (side 
images), shown from the side of the condensate, see Figs. [^, and a. The experimental 


3 



images are carefully compared to the corresponding 2D density maps (Fig. [^) and 3D 
isosurfaces (Fig. [^) resulting from the numerical simulations. 

Given the relatively low condensate temperatures, we model the condensate’s dynamics 
using the Gross-Pitaevskii equation (GPE) [21] 

where t^(r,t) is the condensate’s wavefunction, r the position, t the time, V{r) = + 

ujy'tf -\-uj‘lz^)/2 = m{urr^ +u)zZ^)/‘^ the trapping potential, g = i.'nt^aa/m the strength 
of the inter-atomic interactions and the s-wave scattering length. The normalization 
is = N where V is the condensate volume. The GPE is cast in dimensionless 

form using harmonic oscillator units and solved numerically using the 4th order Runge 
Kutta method with the help of XMDS [H2] • We simulate the vortex decay for a BEG cloud 
with N ^ 1 X 10^ atoms and the same radial and axial harmonic trap frequencies of the 
experiment. 

The top left picture in Eig. shows a typical BEG axial optical depth image of a n = 4 
charged vortex before it decays; on the bottom left, the corresponding number density 
surface. Note the large central circular depletion, representing its core. The top right and 
bottom images show a BEG containing four singly-charged vortices nucleated by the twisted 
decay of the initial quadruply charged vortex. 

Eig. I (a and b) show the computed 3D isodensity surface ([^) on top of its correspond¬ 
ing experimental absorption image (|^), respectively. Note how the vortex lines appear 
intertwined, in the form of helical Kelvin waves along the .s-axis, as predicted [THl EOl 1^ . 
We interpret the central region of Eig. [^, containing a nearly regular density modulation 
along the symmetry axis, as a signature of the presence of Kelvin waves. Eig. (c and 
d) are shown, allowing for a direct comparison of the similar features appearing on the 
experimental (Eig. [^) and the numerical (Eig. [^) column density images. Moreover, we 
present experimental (Fig. [^) and numerical (Fig. [^) optical density profiles, OD(z), di¬ 
rectly extracted from the central (shaded) regions of Figs. and H. along the . 2 -axis. The 
regular oscillations observed near the center on the OD(z) profile demonstrate that the den¬ 
sity modulations are due to the presence of Kelvin waves, clearly visible in the numerical 
3D isosurface (Fig. [^). Likewise, the observation of Kelvin waves reported by Ref. [28] 
was supported by similar evidence. Without Kelvin waves, the central regions of the 2D 
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FIG. 1. Top: Image showing the n — A multi-charged vortex under free fall(left) and the four 
singly-charged vortices resulting from the n = 4 vortex decay (right). Both top images present 
the integrated number density acquired along the BEG z-axis, after f = 14 ms of TOE. Bottom: 
corresponding 3D representation of the top right and left images of the optical densities to empha¬ 
size the vortex core depletions. Typical optical density linear prohles used for data analysis are 
displayed in Eig. 2 of the Supp. Mat. |HH] . 

column density images would be smooth, and the ID OD(z) profiles would not present the 
oscillations shown in Figs, [^e and f), as demostrated in Fig. 3 of the Supp. Mat. [33] . 

The shape of the vortex lines resulting from the decay of multi-charged vortices depends 
on where and when the decay starts. They may appear as straight vortex lines or intertwined, 
as here reported, depending on the perturbation’s symmetry and the local density homo¬ 
geneity. If the perturbation is mostly constant, along z, and the density does not vary much 
in the z-direction, every point on the vortex unwinds at the same rate, and singly-charged 
vortex straight lines are expected to emerge. However, if the density changes significantly 
along z, the unwinding takes place at different times and at different z positions, inducing 
the intertwining, as discussed by [20] . 

The complete time sequence of the twisted vortex decay is revealed by the numerical 
simulations, as shown in the movie [35|. Fig. ([^ presents the three key features (the multi- 
charged vortex, the twisted unwinding, the four singly-charged vortices) observed during 
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FIG. 2. Top: (a) the 3D isosurface resulting from the GP simulation solved for the same exper¬ 
imental conditions, (b) side absorption image showing central density oscillations resulting from 
the twisted unwinding decay. Center: (c) side experimental absorption image, also shown in (b), 
presented in false color alongside its corresponding (d) 2D column density resulting from the GP 
simulation. Bottom: (e) experimental OD z-profile, and the corresponding (f) numerical central 
OD{z) along the condensate 2 ;-axis (averaged over the narrow vertical gap in between the hori¬ 
zontal strips shown). In (e) The smoothed (blue) curve plotted on top the raw experimental OD 
(green) curve after application of a high-frequency filter. 

the intertwined decay at different stages. It is worth comparing this effect, in which helical 
Kelvin waves arise from the interaction of parallel vortices in a confined geometry, with 
the Crow instability [36], which generates Kelvin waves on anti-parallel vortices. Individual 
helical Kelvin waves have been recently observed in superfluid helium following a vortex 
reconnection [37|, whereas in our case, the waves arise from the decay of multiply charged 
vortices, as predicted by [TSl 1^ . 

Numerical experiments suggest that the decay start of the multi-charged vortex can 
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FIG. 3. Numerically computed time sequence of 3D isodensity plots showing the initial n = 4 mnlti- 
charged vortex (a) evolves to the twisted unwinding (b), and finally decays to four singly-charged 
vortices (c). 

be sped up. Imposing random fluctuations 10% of l-^l) to the initial n = 4 vortex 
wave function does not significantly change the decay time scale, probably because the 
symmetry of the initial condition is not completely broken. A small displacement of the 
vortex core axis (~ Uho/S) is more efficient, triggering the onset of the twisted unwinding 
in about 12 ms; a larger displacement (~ Q-ho/2) reduces this time to 10 ms. A few other 
methods were investigated, and the most efficient method found was to gently squeeze the 
harmonic potential in the xy direction by an amount oJxIoJy = 0.9, when preparing the 
initial state in imaginary time, then resetting u)xl<^y = 1 while propagating the GPE in 
real time; this triggers the onset of decay in 5 ms. Experimentally, we found it difficult to 
control perturbations well enough to reproducibly determine the time scale of decay. It was 
observed that, once in inverted bias field configuration, the small instabilities existing on 
the current controller were sufficient to trigger in just 5 ms, which also corresponds to the 
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hold time used for taking the data presented. 



VijOi 

FIG. 4. Turbulent decay of two antiparallel n — 2 vortices. PDFs of turbulent velocity components 
(blue, red, and green symbols respectively for x, y, and z components) at t = 53 ms. Note the 
power law nature of the PDFs at large velocity, unlike the Gaussian fits (gPDFs), which are plotted 
as solid lines. 

As a final remark, it is worth mentioning that the observed twisted decay of multi- 
charged vortices may be exploited to induce quantum turbulence. The existing methods 
used to generate turbulence (rotations j38] , trap oscillations [39] , moving obstacle [iQl l40l - 
1^ . and imprinting staggered vortices [2]) tend to significantly perturb or even fragment 
the condensate. This complicates detecting the vortices position and counting their number, 
which affects the vortex line density estimation. Here we numerically imprint two antiparallel 
doubly quantized (n = 2) vortices in the xy plane. The vortices unwind and twist, moving 
slightly forward due to their self-induced velocity, then the helical waves travelling in opposite 
directions reconnect, generating a turbulent tangle with only moderate density oscillations 
[3H]- We find that the condensate’s velocity field, initially anisotropic {vy/vx ~ 1 and 









Vz/vx ~ 0), becomes almost isotropic in the turbulent stage [vy/vx ~ 0.98 and Vz ~ 0.63), 
displaying velocity PDFs with the typical power-law scaling at high velocity (see Fig. 
observed in superfluid helium [T] and in larger, initially isotropic condensates [2] (whereas 
in ordinary turbulence such PDFs are Gaussian). 

In conclusion, by carefully comparing experimental and numerical results, we have demon¬ 
strated that the decay of multiply charged vortices, topologically imprinted in trapped 
atomic BECs, generate intertwined Kelvin waves which twist and split into singly-charged 
vortices. Numerical experiments suggest that the onset of the decay is sensitive to small 
perturbations. We have also shown that the twisted decay can be used to generate isotropic 
turbulence, relatively free from large scale fluctuations or fragmentation of the condensate, 
hence suitable for comparison with superfluid helium turbulence. Finally, an interesting 
follow-up study would be to assess how the proximity of the BEG edge and the presence of 
the thermal cloud would affect the Kelvin waves [l9]. 
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INGT) and EPSRG. G.D.T and P.E.S.T thank G. Roati for technical support during the 
initial experiment runs. V.S.B thanks E.A.L. Henn for the experimental support. 
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SUPPLEMENTARY MATERIAL 


1. Topological phase imprinting 

The multi-charged vortices are nucleated by adiabatically reducing and inverting the 
initial trap bias field from 5o ~ 0.5 G down to Bq ~ —0.5 G. The resulting magnetic field 
along the weak trap direction, . 2 , axis is reversed during the process. Fig. presents the 
general idea with the main details. During this adiabatic reversal process, the atomic spins 
are forced to follow the rotation taking place in the local field directions. The magnetic field 
reversal rate is slow enough (0.1 G/ms) so that almost all the atomic spins adiabatically 
follow the magnetic field directions, under our typical experimental conditions. Also, since 
the direction of the rotation depends on the spatial position of each atom, the atom acquires 
a different topological phase depending on its position. The spin reversal results in the phase 
winding in the condensate order parameter with dyrh per atom, corresponding to a n = 4 
charge vortex imprinted on the BEG. Finally, this process works well though it induces 
significant loss of atoms (50%) during the zero crossing, due to the Majorana flips, which is 
consistent with previous studies.. 
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FIG. 5. Topological phase imprinting of a multi-charged vortex. The diagram summarizes 
the change in the axial magnetic field profile in the topological phase imprinting procedure used to 
nucleate multi-charged vortices in our BEC. (a) The initial (standard) magnetic field axial profile 
of our QUIC trap, held until the phase imprinting starts, (b): The final magnetic field axial profile 
after the imprinting finishes, about 5 ms later, just before it is switched off for the time-of-flight 
imaging. 


2. Profiles of multi- and singly-charged vortices 

In Fig. typical optical density ID profiles, OD(z), are presented. On the left, in 
Fig.g a), a n = 4 multi-charged vortex with its large central depletion is shown; and, on 
the right, the four singly charged vortices resulting from its decay Fig. [^b) are presented. 
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FIG. 6. Profile of multi- and singly-charged vortices ID optical density profiles corresponding 
to the Fig. 1 of the article, for the multi-charged vortex (left) and the four singly-charged vortices 
(right). 


3. Profile of a vortex without Kelvin waves 


Fig. shows the smooth experimental (b) and numerical (d) ID OZ(z) profiles of the 
central region of the condensate along the z-axis for a vortex without Kelvin waves, corre¬ 
sponding to experimental (a) and numerical (b) false-color absorption images. These smooth 
profiles must be compared to the oscillating profiles for a vortex with Kelvin waves shown 
in Fig. [^c,d,e,f) of the article. 




FIG. 7. Multicharged vortex as seen for short hold times (2.2 ms): (a) experimental ab¬ 
sorption images (false color); and (c) numerical 2D column density, (b) corresponding experimental 
OD 2 ;—profile; and (d) numerical central OD( 2 ;) along the condensate 2 ;—axis (averaged along the 
narrow horizontal strips shown). 
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4. Generation of quantum turbulence 


To exploit the twisted unwinding of the multi-charged vortices as a technique to generate 
vortex tangles, relatively free of large density modulation, we numerically imprinted two 
anti-parallel, doubly quantized vortices, using the parameters listed below, at the initial 
condition, see Figj^a). The vortices unwind, twist, move slightly forward due to the self- 
induced velocity field, and reconnect (Figj^b)), generating a turbulent tangle, see Figgc). 

Initial couditiou: a doubly quantized {n = 2) vortex is imprinted in the xy plane at 
(1.8£r, l-5£r), and a double quantized antivortex is imprinted at {l.Oir, — l.Mr). All param¬ 
eters as in the numerical simulation described in the article. 

Velocity statistics in the turbulent state: the PDFs of the velocity components are 
PDF(ni) ~ n"* {i = x,y,z) with ~ —3.01, ay ~ —3.14, az ~ —3.12. The normalized 
mean values of the Gaussian fits (i.e., mean values gPDF(ni) = ( 27 rcr^)“^A gxp [—(n* — 
divided by their corresponding Gaussian widths) are marked as vertical dashed lines in Fig. 
4. and are respectively given by Vj = —0.014, —0.210, —0.005, with associated widths 
(Tj = 2.3, 2.4,1.5 and /ij = —0.014, —0.210, —0.005 (for i = x,y, z). 
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FIG. 8. Isodensity plots showing the evolution into turbulence of two initial n = 
tiparallel vortices at t = 0 ms (a), t = 21 ms (b) and t = 53 ms (c). 
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